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1. INTRODUCTION 

Today, wireless power transfer (WPT) technology is increasingly attractive for industrial applications 
as well as electric vehicles (EVs) because it simplifies the power supply process and eliminates some of the 
dangers of electrical leakage for users [1]—[3]. Research is toward applying WPT to charge EVs on the move 
because it can reduce the battery's capacity and increase the travel distance [4]—[7]. One of the most important 
criteria in the WPT system is efficiency [8], [9]. Transfer efficiency depends on the compensation circuit, coils, 
working frequency, and load impedance [1]. When the coil, compensation circuit, and working frequency has 
designed, the transfer efficiency only depends on the load impedance. The transfer efficiency is maximum only 
at a load impedance, and it's called the optimal load impedance value [1]. Efficiency drops sharply at other 
load impedance [10]. And, the load impedance changes according to the state of charge of the battery [11]—[13]. 
Therefore, impedance matching control is necessary to improve efficiency. 

In addition, the optimal load impedance depends on the coupling coefficient. In dynamic wireless 
charging systems, the coupling coefficient varies with the position of EVs [14]—[16]. Therefore, the problem 
of estimating the coupling coefficient when the vehicle is moving needs to be done. Many studies have been 
carried out to estimate the coupling coefficient in the static charging system [17], [18], as well as in the 
dynamic charging system [19]—[22]. However, these methods only work well under ideal conditions. Under 
actual operating conditions, the characteristics of the passive elements may be changing, measurements may 
be affected by noise signal. These conditions make the above estimation methods inaccurate. 
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The Kalman filter is used to estimate a state variable that is not measurable, or measurable but inaccurate 
due to the impact of measured noise with high accuracy [23], [24]. Therefore, to overcome the above 
disadvantages, this paper has proposed a coupling coefficient observer based on the Kalman filter. Firstly, the 
designed structure of the dynamic charging wireless system is presented. Secondly, the improving transfer 
efficiency solution is analyzed. Thirdly, a coupling coefficient observer base on the Kalman filter is built. Finally, 
on the secondary side, an optimal load impedance controller is realized using an active rectifier. The results show 
that the coupling coefficient is estimated with high accuracy, and the efficiency system is improved. 


2. EFFICIENCY-IMPROVING CONTROL BASE 
2.1. System structure 

Figure 1 shows the proposed dynamic wireless charging system for EVs. The primary side is the 
modular design, each transmitting module includes a high-frequency inverter that powers the three transmitting 
coils via LCC compensation circuits. The coils of the transmitting modules are arranged under the roadway to 
form a dynamic charging lane for EVs. On the secondary side, the receiving coil is mounted under the EVs to 
receive energy from the dynamic charging lane. An active rectifier is used for optimal impedance control to 
maximize transfer efficiency. And then, the energy is passed through the battery energy management controller 
to charge the battery. This paper only considers the problem of wireless power transfer on a transmitting 
module, not considering issues such as the switching between transmitting modules and battery energy 
management control. 
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Figure 1. The proposed dynamic wireless charging system 


2.2. Solutions to improve transfer efficiency 

The LCC compensation circuit has many advantages and is suitable for dynamic charging 
systems [25]. The double-sided LCC compensation circuit was selected for this system design. Figure 2 shows 
the equivalent circuit. Rz is the equivalent load impedance from the active rectifier to the battery. The 
electromagnetic coupling between the coils is represented by the induced voltages. The induced voltage in each 
coil depends on the value of the mutual inductance and the current in the other coil. Each transmitting coil 
connects electromagnetically to the other two transmitting coils, and it is replaced by a jæM;Iz; induced voltage. 
Each transmitting coil connects electromagnetically to the receiving coil and it is replaced by a jæMi-Lz, induced 
voltage. The receiving coil connects electromagnetically to the three transmitting coils and it is replaced by 
joMJ induced voltage [25]. The transfer efficiency can be calculated as (1) [25]. 


R,I? R 
L'Lfr _ L (1) 


1) = S072 2 2 2 z= 212 
Ri, pet Rrliyp + Ralf +R212+R31f3 2_Rr Sen (14 6 ) 3o“ Ltr 1 
L wL, k2.Q;Qr k2.Q;Qr Rr KẸQiQr 


Where, Q; = œL/R;, Q, = wL/R, are the quality factors of the transmitting and receiving coils, respectively; k, 
is the coupling coefficient. The (1) shows that with a designed system working at the resonant frequency, the 
transfer efficiency depends on the equivalent impedance load. Solving in (2) obtains the maximum transfer 
efficiency as shown in (3). 


2 
ŽL =Oand 5 <0 (2) 
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The maximum transfer efficiency is determined by (3). 
= kK? QiQr 8) 


Nmax = 2 
(va f+ržoior) 


Is achieved at (4). 


wL? 
— Pein = 
Ropt ~ Ry 3+Kk7QiQr (4) 


The (3) and (4) show that when the equivalent load impedance (Rzo);) is the optimal load, transfer 
efficiency is maximum. According to (4), the optimal load depends on the coupling coefficient. However, in a 
dynamic wireless charging system, the coupling coefficient depends on the EV's position. Therefore, an 
estimation of the coupling coefficient value is necessary to improve transfer efficiency. 

To estimate the coupling coefficient, perform an equivalent circuit analysis of the system at resonance 
condition in Figure 2. Study [22] has presented a method to estimate the coupling coefficient using the (5). 


= a Rr 
r Uap Jlilr lol yr Uap + wLprlifr 7 


According to (5), if the r.m.s. of the voltage (Uap) and current (/7/,) are measured, the coupling coefficient can 
be estimated. However, (5) is established under ideal conditions, if the parameters of the passive elements 
change, the accuracy of the estimation method decreases. 


joMAi 


Figure 2. Equivalent circuit 


2.3. Building the coupling coefficient observer based on the Kalman filter 

The target of controlling the load impedance according to the optimal load is only possible when the 
coupling coefficient is known. However, the conventional estimation method only works well under ideal 
conditions. Under actual operating conditions, the characteristics of the passive elements may change, which 
reduces the accuracy of the estimation method. This section proposes a method to observe the coupling 
coefficient based on the Kalman filter, which overcomes the above disadvantages. 


2.3.1. System modeling 
The induced voltage on the receiving coil is performed as follows. 


E = oMyly, ©) 


Where M, is mutual induction of transmitting coils with receiving coil (M, = k,../L;L,). 
The resonant current in the transmitting coils is the same [25]. 


Ii = —joCp Uap (7) 


Where Uag is the r.m.s of the output voltage inverter at the primary side. In (7) shows that the resonant current 
on the transmitting coil is independent of the parameters on the primary side. When the transmitting power was 
kept fixed, the resonant current is constant. Therefore, the mutual induction (M,) can estimate by estimating 
the value of the induced voltage (E) as in the (8). 
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a. (8) 


wlyi 


On the secondary side, the Kirchhoff equations are written for the equivalent circuit in Figure 2 as 
follows: 


di P 
L, = + Ucr + Ucrr = E sin(wt) 
ducr _ i 
r = REP 
ae ©) 
Cfr ae = lir — lifr 
diz f 3 
Ler a + typ Ry = Ucfr 


Using the first harmonic approximation (FHA) [26], convert current and voltage values into sine and cosine 
components: 


x = Xsin(wt + p) = x, cos(wt) + x, sin(wt) (10) 
Derivative of sine and cosine components: 


dx dx¢ dxs 
Fs (= + wxs) cos(wt) + (= + wx) sin(wt) (11) 


Rewrite the (9) in terms of cosine and sine components: 


di 7 
Ly ( ae + wirs) + Ucr,c + Ucfric = 0; 


dizys . - 
at dt = Wlirnc + Ucr,s + Ucfr,s =E 
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Cal dt + wucr,s) cs lLr,c) 
ducr,s ne 
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Lfr ( dt + witprs) + ingreRi = Ucfr,ci 
diLfr,s . is 
Lfr ( dt + wirpre) + iLfrs RL = Ucfr,s 
State variables are set as shown in Table 1. The state model represents the kinematics of the system: 
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Consider E to be the 9" state variable (xo) and to be fixed relative to the rest of the system's state 
variables at the time of consideration. 


4, ===0 (14) 


dt 


Combining (13) and (14), the system model is as follows: 
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x = Ax (15) 
Where: 
0 -w -0 0 -6 0 0 0 0 i 
w 0 0 -6% 0 -6 0 0 6 => 
6, 0 0 -w 0 0 0 0 0 A 
0 & ow 0 0 0 0 0 0 =. 
A=|@6; 0 0 0 0 -w -0 0 O|. with D (16) 
0 6 0 0 w 0 0 <6. 0 oars 
0 0 0 0 8% 0 -R,6, -w 0 4 
0 0 0 0 0 & w RO, 0 8, = fee 
0 0 0 0 0 0 0 0 0 


The parameter Rz changes during charging, so the kinetic model is nonlinear. Implementing the 
discretization algorithm for each control cycle, the system model will have the form: 


Xk41 = Xg; with Ọ = e's4 (17) 


Where 7; is the calculation period of the Kalman filter. 


Table 1. State variable name table 


Parameter Value 
Variable line lirs Ucr,c Ucr,s Ucfric Ucfrs lLfr,c lLfr,s 
State EA Xz X3 X4 Xs X6 X7 Xg 


2.3.2. Apply the Kalman filter to observe the coupling coefficient 
The measurand is the r.m.s. of the rectifier input current it is expressed in state variables as follows. 


Y = lişr = hx) = x7 + xg (18) 
The Kalman filter consists of two main calculation stages: 


- Stage 1: Prediction. Based on the mathematical model, the Kalman filter calculates the values of the next 
state variables (xz+1) and the deviation of those calculations from reality (Px+1). 


ia = ỌkXk (19) 
f — T 
Pear = PrP Pi + Q (20) 
- Stage 2: Adjustment. The measured signal had sent to the microcontroller (with measurement noise). At 


this stage, the Kalman gain is updated, and the state variable vector (xk) and its deviation (Pk) are correct. 
The Kalman gain, the state variable vector, and the covariance matrix had corrected according to (21)-(23). 


Ky = PI HE (H,P{ HE +R) (21) 

Xe = Xf + Ke(Ye — hCxf)) (22) 

Py = (I — Ky Hy) Pf (23) 
Were, Hk = gn) 2X7 + 2xg. Hy is the Jacobian matrix of h(x), and Hx is updated after each calculation 


dx 
cycle. Q and R are covariance matrices of the measurement and process noise. 


The Kalman filter algorithm diagram is shown in Figure 3. The current in front of the active rectifier 
(Ij) is measured and fed into (22). Kalman filter is based on the state model and the values of state variables 
in the past and calculates the vector of state variables at present (xz). The optimal load impedance is estimated 
based on the coupling coefficient estimation results (4). Then this value is discretized to calculate the matrix A 
(16). The state variable is calculated according to (19). the Kalman gain is updated, and the state variable vector 
(xx) and its deviation (Px) are correct. 
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Figure 3. Kalman filter algorithm 


2.4. Analysis of impedance conversion method using an active rectifier 

Usually, on the secondary side, a diode bridge rectifier and DC/DC converter are used to impedance 
conversion. This paper proposes to replace it with an active rectifier. This solution reduces one inductance and 
one filter capacitor. As a result, the size, weight, loss, and cost of the receiver are reduced. Moreover, this paper 
uses the phase synchronization method combined with a symmetrical phase shift to convert load impedance. 
The input voltage and current of the active rectifier are phase-synchronized so that the equivalent impedance 
is purely resistive. Then perform a phase shift to control the load impedance value. 

For simplicity of analysis, consider the MOSFET is ideal, ignoring the losses on the filter capacitor 
Co and the capacitor Co has a large enough value for the output voltage Uze to be constant. Rze is the equivalent 
load impedance from the battery energy management controller to the battery. Figure 4(a) is a waveform in an 
active rectifier. Where, u'a» and i’, are the fundamental harmonic component of the input voltage and current 
of the active rectifier, respectively. Vo; to Vo4 are signal turns on MOSFET from Q; to Q4. The MOSFET turn- 
on signal has synchronized with the moment when the input current (i’,») is zero. In each cycle, there are four 
switching states and six phases, as shown in Figure 4(b). 
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Figure 4. Active rectifier (a) control signal and (b) on/off state of MOSFETs 


State A, corresponding to stages (1) and (3): turn on Q2 and Q4. Current ia flows through Q2 and Q4, 
and energy is not transferred to the load. Therefore, voltage ua» and current Ize are zero. Energy has accumulated 
in the capacitor Co fed to the load. State B, corresponding to stages (1) and (3): turn on Q; and Q2. Current iab 
flows through Q; and Q72, and energy is transferred to the load. The analysis is similar for states C and D. Thus, 
in a cycle, uap has the values in the (24). 
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0:0 < wt < ß/2 

Ure: B/2 < wt <r -—pfp/2 

0O: —p/2< wt <a+B/2 
—Ue: T+ B/2< wt < 2r —fß/2 


tis = (24) 


Where £ is the phase-shift angle; Uze is the output DC voltage of the active rectifier. Based on the fundamental 
harmonic analysis method, the r.m.s. of uap is calculated as 25. 


Us = u,-cos? (25) 


w 
Based on the power balance condition, Rz can be calculated as in (26). 
8 B 
R, = zz Riecos* a (26) 


In (26) shows that the equivalent load impedance (Rz) can be adjusted by changing the phase shift angle. 

Figure 5 shows the block diagram of the optimal impedance controller structure. The r.m.s. of the 
active rectifier input voltage and current are measured to calculate the equivalent load (Rz) and estimate the 
optimal load impedance (Rzp:) according to the Kalman filter algorithm. The error value is fed to the PI 
controller. The PI unit calculates the phase shift angle for the active rectifier. In addition, the PLL is used to 
determine when the current i,, is zero, and to generate a synchronous signal that is fed to a phase shift signal 
modulator that drives the MOSFETs. 
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eee modulation Q3 
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Figure 5. Optimal impedance controller block diagram 


3. RESULTS AND DISCUSSION 
3.1. Kalman filter verification 

A system simulation model is built on MATLAB/Simulink software to verify the accuracy of the 
Kalman filter. The system parameters are built according to Table 2. The value of the mutual inductance 
between the transmitting and receiving coils is approximated by the (27). 


M = Map + 20%. Mypsin(2nfyt) (27) 


Where M,, = 16.4443 uH, with a vehicle speed of 40 km/h and a transmitting coil length of 40 cm, M changes 
with frequency f= 27 Hz. The mutual inductance characteristic is shown in Figure 6(a). The simulation model 
considers the impact of measurement noise as in Figure 6(b). The noise source has a frequency of 85 kHz and 
an amplitude of 10%. The actual and estimated mutual inductance values for several cases are in Figure 7. In 
case, the system parameters are equal to the original design parameters. The simulation results in Figure 7(a) 
show that the estimated inductance value follows the real, with an error of less than 5%. In case, system 
parameters change from the original design parameters. The parameter's value varies according to Table 3. 
Figure 7(b) shows the estimated results in the three cases as follows. 

In case 1, regardless of the impact of measurement noise, the mutual inductance (M) is estimated 
according to (5). The results show that the estimated error is 14.5%. In case 2, considering the impact of 
measurement noise, M is estimated according to (5). The results show that the value of M received remains the 
same as the amplitude of the measured noise, and the estimated error is over 14%. In case 3, considering the 
impact of measurement noise, M is estimated according to the Kalman filter algorithm. The results show that 
the value of M follows the actual with an estimated error of less than 5%. Thus, the results of estimating M by 
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the Kalman filter algorithm have high accuracy in the case of measurement noise impact and the system's 
parameters change. 


Table 2. System parameters 


Parameter Value Parameter Value Parameter Value 
Po 1.5 kW L, 120 pH C2 123.2 nF 
Upc 310 V R, 0.13 Q C3 95 nF 
Few 85 kHz Li 52.6 yH Ly 28.9 uH 
Li 102 H Ch 66.5 nF Ch 120.9 nF 
Ri 0.13 Q Gi 93.7 nF C, 38.5 nF 


7 7 
[= Noise 

Mutual j 5.5 |—Real value |] 

inductance ‘ i T 


Current 
If ( 


i L i | i i | i L 
0.018 0.02 0.022 0.024 0.026 0.028 0.03 0.032 0.034 
Time (s) 


(b) 


Figure 6. The characteristic (a) mutual inductance approximation and (b) measurement noise 


Table 3. System parameters change from the original design parameters 


Parameter Li Li Ci C; Ri Ly R, 
Rate -5% (down) -10% (down) -10% (down) -6% (down) 15% (up) -7% (down) 8% (up) 
x10 x105 
i 7 
2+ — — J E i J 
g F F me F Ta g ’ F 
s 187 {f N A AS] 4 = L8p = J 
S = 7 = H 1A 3 4 
i 1. l! z NI ERK ~~ 
314} = Í 5 EEP s ~ 
f=] ap k=] " 
E a i 
2127 J 2 1.2} - 
Š = M_real (preset) 
14 M real (preset) ||| |—M_estimation- in case 1 
a Hp mm ~~ |-—M_estimation- in case 2 
08 ji ji i 1 == oR | | | |—M_estimation- in case 3 
0 0.01 0.02 0.03 0.04 0.05 0.06 “0 0.01 0.02 0.03 0.04 0.05 0.06 
Time(s) Time(s) 
(a) (b) 


Figure 7. Mutual inductance estimation results (a) original design parameters (b) system parameters change 


3.2. Experimental results 

An experimental model of dynamic wireless charging with a power of 1.5 kW, transfer distance of 
150 mm, and working frequency of 85 kHz was built in the laboratory to verify the proposed method and is 
shown in Figure 8. The C3M0280090D MOSFET SICs are used to increase the converter efficiency. 
Polypropylene film capacitors are used as compensating capacitors because of their low loss and high- 
frequency current tolerance. Ferrite bars are used to increase the electromagnetic connectivity between the 
transmitting coils and receiving coil. 

Figure 9(a) shows the measurement results of the input and output signal of the phase synchronous 
circuit. The results show that the output pulse of the PLL circuit has accurately captured the phase angle of the 
current. The pulse is high when the current cuts through zero from negative to positive and is low when the 
current cuts through zero from positive to negative. Active rectifier input voltage/current are shown in 
Figure 9(b) which corresponds to the case of phase shift angle of 30 degrees. The results show the active 
rectifier's ability to control phase shift. 
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When the equivalent load impedance varies between 80 Q+130 Q, the load impedance response is 
shown in Figure 10(a). The optimal impedance (RLop) is obtained from the results of the Kalman filter 
algorithm. Although the EV is fixed position, the estimated optimal load impedance changes slightly. It is due 
to the effect of current and voltage measurement noise. The fluctuation amplitude of the optimal impedance is 
estimated to be no more than 1%. The results show that the equivalent impedance (Rz) follows the optimal 
impedance (R_z opr) and the static error is less than 1%. 
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Figure 8. Experimental model 


OFreqency 64.95kii2 Max 7.600 Ørin -7.400 
@Max 4.720 @Min -488nV @Pk-Pk 5.120 


82 p---- aS) ene lea aaa cat eas ania seca 1 


oo 
Oo 


9 


78 e $- a o--- ee or > -È i -$ 


1 1 
1 1 1 
rererere pose 
1 1 1 


System efficiency (%) 
= 
a 


52.5 p~ ——RL_Equivalent impedance (response) g ' ' ' ' 
kipi 72 Peco Arcencie sco bees seee@ eee Eff, Case2 ===! 
= RL.opt_Optimal impedance ! i i @— Eff Case3 i 
52 1 1 1 i i i i i i 70 ! i H i i i 
O 10 20 30 40 50 60 70 80 90 100 0 100 200 300 400 500 600 700 800 
Sample cycle (kTs), Ts=1 1.7us Position's EV (mm) 
(a) (b) 


Figure 10. Characteristic (a) load impedance response and (b) system efficiency 


When EV moves along a dynamic wireless charging lane, the system efficiency experimental results 
from the primary inverter input to the equivalent load (Rze) are shown in Figure 10(b). System efficiency 
includes inverter efficiency, transfer efficiency, and active rectifier efficiency. 
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- Case 1: system parameters are equal to designed parameters as in Table 2. When the optimal impedance is 
estimated using the Kalman filter algorithm, the average system efficiency is 80.76%. 

- Cases 2 and 3: system parameters are changed as shown in Table 3. When the optimal impedance is 
estimated using the Kalman algorithm, the average system efficiency is 77.94%. When the optimal 
impedance is estimated using (5), the average system efficiency is 74.71%. In this case, the system 
efficiency is improved by 3.2% when using the Kalman filter algorithm. 


4. CONCLUSION 

The paper proposes to observe the coupling coefficient based on the Kalman filter for dynamic 
wireless charging systems. The improving transfer efficiency solution is analyzed. The optimal impedance has 
been estimated through the coupling coefficient estimation, with a static error of less than 1%. Then, the phase 
synchronous and phase shift methods are used to control the optimal impedance. The results show that the 
efficiency is improved, and the experimental system efficiency is 80.76%. 
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